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Recursion relations for Hylleraas three-electron integral are obtained in a closed form by using 
integration by parts identities. Numerically fast and well stable algorithm for the calculation of 

^ | the integral with high powers of inter-electronic coordinates is presented. 
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I. INTRODUCTION 

The explicitely correlated Hylleraas basis set Q is one of the most efficient representation 
of a few-electron wave function. It has been applied mostly to helium and lithium atoms 
3]. While the general analytic formula for two-electron Hylleraas integrals is quite simple 
4], the three-electron case is much more complicated. Fromm and Hill in j^] were able to 
derive a closed form formula for the generic integral 

J d 3 n J d 3 r 2 J d 3 r 3 e-^iri-iwra-toars-uiras-iiarsi-iisria r ~\ r ~l r ~l r ~l r ~l r ~l ^ 

which consists of the sum of many multi-valued diloearithmic functions. This formula was 

n 

later simplified by Harris [6] to a more symmetric form, in which many spurious singularities 



have been eliminated. However, only few preliminary results have been reported so far with 
these tati o„el yinM , 03. The more convenient basis set for applications to 
3- or 4-electron atoms happen to consist of functions with all itj equal to 0, and include 
polynomials in and additionally in r^, namely the 3-electron Hylleraas functions 



r , _ „-w 1 ri-102 r 2 -w 3 r s ni-1 n 2 -l „n 3 -l n 4 -l r n 5 -l „n 6 -l 
{n} — e ' 23 ' 31 ' 12 ' 1 '2 '3 



(2) 



A compact expression for the integral with all = and with n\ = 2, n 2 = 0, n 3 = 2, = 
= uq = 1 was derived in jjj. No attempts have been reported so far to derive analytic 
formulas for larger values of rij. They can in principle be obtained from the general formula 
by Fromm and Hill in Ref. p by differentiation with respect to Ui and Wi. In practice 
that is quite complicated for large powers of r^. For testing results obtained in this work 
we have written a Mathematica code for derivation of Hylleraas type of integrals using the 
Fromm-Hill formula, and obtained results for all rii such that n\ + n 2 + ^3 + ^4 +72.5 + n e < 10. 
Despite the existence of an analytical result for the generic integral (0), the most advanced 
results for lithium have been obtained so far by Yan and Drake in series of papers |3||, where 
Hylleraas integrals where calculated by series expansion in angular momenta. 

For helium atom recursion relations J^J are available, which allow to express Hylleraas 
integrals with n\ + n 2 + ^3 = N [see Eq. ()A5)I ] as a linear combination of integrals from the 
lower shell N — 1. In this work we derive a similar complete set of recursion relations for 
the three-electron Hylleraas integrals 

/Kn 2 ,n 3 ;n 4 ,n 5 ,n 6 ) = [ ^ [ ^ [ ^ — — 

J 4lT J 4 7T J 4lT 

r ni— 1 „n2-l -713-1 r nA—\ r ns— 1 -ne-l (o\ 
'23 '31 '12 '1 '2 '3 • \°) 



These recursions allow to eliminate the principal drawback of Hylleraas basis sets, namely 
the complexity of various integrals for nonrelativistic matrix elements for the lithium atom. 
Moreover they can be applied also to more than three electron atoms with the condition 
that the basis set contains at most one odd power of r^. In the next section we present our 
main result, a complete set of recursion relations for the three-electron Hylleraas integrals 
defined in Eq. (J3J). In Section III we analyze boundary terms and two-electron recursions, 
and in Section IV we sketch the derivation of the main result. Finally in Section V we 
analyze numerical stability and present numerical results for some selected integrals. 



II. THREE-ELECTRON RECURSION RELATIONS 



The recursion relations are conveniently divided into two sets. The first set of three 
recursions increases by 2 any of the first three indices ni,n 2 , n 3 at n 4 = = = 0, starting 
from the shell n\ + n 2 + n 3 = (as initial conditions, the values with n« = 0,1,2 = 1,2,3 
and 77.4 = riz, = n§ = are therefore needed); with the notation of Appendix A, they read 

1+^3 



/(n 1 ,n 2) n 3 + 2; 0,0,0) 



n 2 (n 2 - 1) 



/(n 1) n 2 -2,n 3 + 2;0,0,0) 



^3 + 1 

+(m + 2 n 2 + n 3 + 2) f(m, n 2 , n 3 ; 0, 0, 0) 

H TT n 2, n 3 + 2; *, 0, 0) 

n 3 + 1 

n 2 (n 2 - 1) 



/(m + 2,n 2 -2,n 3 ;0,0,0) 



ni + 1 

+ — r f(ni + 2, n 2 , n 3 ; 0, 0, ★) 

Tl\ + 1 



'n 2 



n 3 + l 
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ni + 1 



/(m,*, 7i 3 + 2; 0,0,0) 
/(n 1 +2,*,n 3 ;0,0,0) 



+ 



Ul ( ^ 1 1) /(n 1 -2,n 2 ,n 3 + 2;0,0,0) 
n 3 + 1 

+ (2 m + n 2 + n 3 + 2) /(m, n 2) n 3 ; 0, 0, 0) 

+ — l —r f(ni, n 2 , n 3 + 2; 0, *, 0) + 
^3 + 1 

ni (ni - 1) 



n 2 + l 



/(n 1 -2,n 2 + 2,n 3 ;0,0,0) + 
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/(n 1 ,n 2 + 2,7i 3 ;0,0,*) 



n 2 + l 



'"' /(*,n 2 ,n 3 + 2;0,0,0) 



n 3 + l 



" ni T /K^2 + 2,n 3 ;0,0,0) 



^3 f ra 3 (n 3 - 1) 



«2 + 

/(m, 112 + 2,713-2; 0,0,0) 

n 2 + 1 

+(ni + n 2 + 2 n 3 + 2) f{n 1 ,n 2 , n 3 ; 0, 0, 0) 

H 7T f(ni, 2 + n 2 , n 3 ; *, 0, 0) 

n 2 + 1 

+ n3 (n3 ,~ 1) /(ni + 2, n 2 , n 3 - 2; 0, 0, 0) 

711 + 1 

+ — r /Oi + 2, n 2 ,n 3 ; 0,*0) 

Tl\ + 1 

/(m,n 2 + 2,*;0,0,0) 



n 2 + 1 

"" ; /(m + 2, Tia,*; 0,0,0) 



+ ' ---- — /jj ( 4 ) 

Formulas for /(ni, n 2 + 2, n 3 , 0, 0, 0) and f{n\ + 2, n 2 , n 3 , 0, 0, 0) can be obtained from Eqs. 
© and (JU) by interchanging (r 2 ,r 3 ) or (ri,r 3 ) respectively. 

The next three recursions increase by 1 any of 714, 715, ri6 of a given subshell n 4 + + 7i6 
for arbitrary m, n 2 and ri 3 

f(n 1 ,n 2 ,n 3 ;n 4 ,n 5 ,n 6 + 1) = { 

iui u> 2 u> 3 

(ni - 1) 7i! n 4 f(m - 2, n 2 , n 3 ; n 4 - 1, n 5 + 1, n 6 ) 

+ (n 2 - 1) n 2 n 5 /(m, n 2 - 2, n 3 ; n 4 + 1, n 5 - 1, n 6 ) 

-(n 3 - 1) n 3 n 4 /(m, n 2 , n 3 - 2; n 4 - 1, n 5 + 1, n 6 ) 

-(ra 3 - 1) n 3 n 5 /(m, n 2 , n 3 - 2; n 4 + 1, n 5 - 1, n 6 ) 

+n 4 n 5 (m + n 2 - n 3 + n 6 + 1) /(m, n 2 , n 3 ; n 4 - l,n 5 - l,n 6 ) 

-(m - 1) Tii /(m - 2, n 2 , n 3 ; n 4 , n 5 + 1, n 6 ) w x 

+ (n 3 - 1) n 3 f(m, n 2 , n 3 - 2; n 4 , n 5 + 1, n & ) w\ 

-n 5 (ni + n 2 - n 3 + n 6 + 1) /(m, n 2 , n 3 ; n 4 , n 5 - 1, n 6 ) »i 

-(n 2 - 1) n 2 /(tii, ™ 2 - 2, n 3 ; n 4 + 1, n 5 , n 6 ) w 2 

+ (n 3 - 1) n 3 /(m, n 2 , n 3 - 2; n 4 + 1, n 5 , n 6 ) u> 2 

-n 4 (ni + n 2 - n 3 + n 6 + 1) /(ni, n 2 , n 3 ; n 4 - 1, n 5 , n 6 ) w 2 



-n 4 n 5 /(ni, n 2 , n 3 ; n 4 - 1, n 5 - 1, n 6 + 1) w 3 

+(m + n 2 -n 3 + n 6 + l) f(n h n 2 , n 3 ; n 4 , n 5 , n 6 ) w 2 

+n 4 /(m, n 2 , n 3 ; n 4 - 1, n 5) n e + 1) w 2 w 3 

+n b f(m, n 2 , n 3 ; n 4 , n 5 - 1, n 6 + 1) w 3 wi 

S ni n 4 /(★, n 2 , n 3 ; n 4 - 1, n 5 + 1, n 6 ) 

+<5 ni /(*, n 2 , n 3 ; 7i 4 , n 5 + 1, n 6 ) u>i 

-^n 2 n 5 /(ni,*, n 3 ; n 4 + 1, n 5 - 1, n 6 ) 

+5 na /(ni ! *,7i3;ra 4 + l,n 5 ,n 6 )w 2 

+^n 3 ra 4 /(ni, n 2 , ★; n 4 - 1, n 5 + 1, n 6 ) 

+^n 3 n 5 /(ni, n 2 , *; n 4 + 1, n 5 - 1, n 6 ) 

-*na /( n i> n 2, *; n 4 , n 5 + 1, n 6 ) t«i 

-^n 3 /(ni, n 2 , *; n 4 + 1, n 5 , n 6 ) w 2 } . (5) 

Two other recursions for f(ni, n 2 , n 3 , n 4 , 715 + l,n 6 ) and /(ni, n 2 , n 3 , n 4 + l,ra 5 ,7i 6 ) are ob- 
tained from Eqs. (JSJ) and (0) by interchanging arguments of /-function, using the identities 

/(ni, n 2 , n 3 ; n 4 , n 5 , n 6 ; w x , w 2 , w 3 ) = 
f(n 2 , ni, n 3 \ n 5 , n 4 , n 6 ; u> 2 , u>i, u> 3 ) = 
/(n 3 , n 2 , Tii; n 6 , n 5 , n 4 ; u> 3 , tu 2 , u>i) = 

/(ni, n 3 , n 2 ; n 4 , n 6 , n 5 ; Wi, iu 3 , iu 2 ) , (6) 
where any of rii can become a *. 

III. BOUNDARY TERMS 

These recursion relations in Eqs. and (J^J) involve boundary terms. The first boundary 
term was derived in Ref. j9|, the 7 remaining boundary terms to start the first set of 
recursions are obtained from the generic Fromm-Hill formula p or by direct integration 

/(0, 0,0; 0,0,0) = — — { 

2 Wi w 2 w 3 I 

U-!*-) h,(i + -^-) + u 2 (—^-) +uJi 

\wi + u' 2 / v wi + w 2 / v wi + w 2 y v u>i + w 2 / 



5 



+ ln(^- 

\W 3 + Wi 
Wi 



+ ln 



7(1, 0,0; 0,0,0 
/(0, 1,0; 0,0,0 
/(O, 0,1; 0,0,0 
/(l, 1,0; 0,0,0 
7(1,0,1; 0,0,0 
7(0, 1,1; 0,0,0 
7(1, 1,1; 0,0,0 



W 2 + W 3 
1 



mil + 
lnf 1 + 



w 2 



Wi + w 3 
Wi 



w\ w 3 



wj w 2 



In 



In 



In 



W\ + w 2 , 
Wi (wi + w 2 + w 3 ) 
(wi + w 2 )(w 1 + w 3 ) 
W 2 (Wi + w 2 + w 3 ) 
(w 2 + w 3 )(w 2 + Wi) 
w 3 (wi + w 2 + w 3 ) 
(w 3 + wi)(w 3 + w 2 ) 



+ Li 2 
+ Li 2 



w 2 



Wi + w 3 

Wi 

w 2 + w 3 



+ Li 2 ( 1 
+ Li 2 ( 1 - 



w 2 



Wi + w 3 
w 2 + w 3 



Wi w 2 (wi + w 2 ) w 3 ' 

1 



2 ■ 



Wi w 3 (wi + w 3 ) w 2 

1 

w 2 w 3 (w 2 + w 3 ) wl ' 

1 



2 2 2 ' 
Wl Wo Wn 



(7) 



J l ^2 ^3 

where Li 2 (x) is the dilogarithmic function. Besides the initial terms one needs also as 
boundary terms the following two-electron integrals 



7(*, n 2 ,n 3 ; 


n 4 , 


n 5 ,n 6 ) 


= T(n 5 + n Q - 


1, n 4 , n 3 + n 2 - 1; w 2 + u> 3 , w h 0) , 


7(ni,*,n 3 ; 


n 4 , 


n 5 ,n 6 ) 


= T(n 4 + n 6 - 


1, n 5 , m + n 3 - 1; w a + w 3 , w 2 , 0) , 


7(m,n 2 ,*; 


n 4 , 


^5,^ 6 ) 


= T(n 4 + n 5 - 


1, n 6 , ni + n 2 - 1; wi + w 2 , w 3 , 0) , 


f(ni,n 2 




★,0,0) 


= T(n 3 - l,n 2 


- l,ni; w 2 ,w 3 ,0) , 


f(ni,n 2 




0,*,0) 


= r(ni - l,n 3 


- l,n 2 ;u; 3 ,wi,0) , 


f(ni,n 2 


% 


0,0,*) 


= r(n 2 - i,m 


- l,n 3 ;wi,w 2 ,0) , 



(8) 

where T is defined in Eq. (IA5I) . Recursion relations for T have already been worked out in 

nn 

the literature and can be found for example in |4|, HQ] • For completeness we include below 
generic formulas for two cases: nonnegative ni,n 2 ,n 3 , 

d \ ni / d \" 2 / d \" 3 1 



T(ni,n 2 ,n 3 ;a u a 2 ,a 3 ) 



dai 



da 5 



da 3 7 (ai + a 2 )(a 2 + a 3 )(a 3 + ai) 



and the second case: rii = —1, and nonnegative n 2 , n 3 

d \ n2 f d \" 3 m(o:i + a 2 ) - \n{a x + a 3 ) 



r(-l, n 2 , n 3 ; a x , a 2 , a 3 



da-: 



da-: 



(a 2 - a 3 )(a 2 + a 3 ) 



(9) 



(10) 
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The complete set of recursion relations have been checked by comparison with analytic 
expression obtained by differentiation of Fromm and Hill formula [5] with respect to «j and 
Wi parameters (in our notation). 



IV. DERIVATION OF THREE-ELECTRON RECURSION FORMULAS 



We use the method of the integration by parts identities [11], which is by now standard 
in the analytical calculation of Feynman diagrams. In our case, it amounts to consider the 
following 9 identities in the momentum representation of the integral G, [see Eq. (jA6j) ] 

d_ 

dki 

(kj + uir m * (ki + uir^(ki 2 + wir^ (** + ^ + ^)- m6 l , (n) 



o = id(i,j) = j d 3 h j d 3 k 2 j d s k 3 -i (ki 



which are trivially valid (the integral of the derivative of a function vanishing at infinity 
vanishes) . 

These identities group naturally into three sets. The first set consists of id(l, 1), id(2, 1), 
and id(3, 1). Other sets are obtained by changing the second argument from 1 into 2 or 3. 
The reduction of the scalar products from numerator leads to the following identities of the 
first set 

id(l, 1) = — m 6 G(mi — 1, m 6 + 1) — m 5 G(m 1 — 1, m 5 + 1) + m 6 G{m 2 — 1, m 6 + 1) 
+m 5 G(m 3 -l,m 5 + l) + (3-2m 1 -m 5 - m 6 ) G() + m 6 G(m 6 + 1) u\ 
+m5 G(m^ + l)?ii + 2mi G(mi + 1) u\ — ttiq G(me + 1) u\ — m.5 G(rri5 + 1) u\ 
+m 5 G(m 5 + 1) w\ + m 6 G(m 6 + 1) w\ , 



id(2, 1) = — m 6 G{mi — 1, m 6 + 1) — m 5 G{mi — 1, m 5 + 1) + m 6 G(m 2 — 1, m 6 + 1) 
+m 5 G(m 3 — 1, m 5 + 1) — m 5 G(m 4 — 1, m 5 + 1) + (m 6 - mi) G() 
+777.5 G(m5 + 1, m 6 — 1) — mi G(mi + 1, m 2 — 1) + mi G(mi + 1, m 6 — 1) 
+m 6 G(m 6 + 1) -u^ + m 5 G(m 5 + 1) u\ + mi G(rrii + 1) u\ — m 6 G(m 6 + 1) u\ 
+mi G(mi + 1) u\ — m 5 G(m 5 + 1) -Ug + m 5 G(m 5 + 1) — m 6 G(m e + 1) u>3 
— m 5 G(m 5 + l)wl — mi G(mi + 1) wi , 
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id(3, 1) = — m 6 G{mi — 1, m 6 + 1) — m 5 G(mi — 1, m 5 + 1) + m 6 G(m 2 — 1, m 6 + 1) 
+m 5 G(m 3 - 1, m 5 + 1) — m 6 G(m 4 - 1, m 6 + 1) + m 6 G(-m 5 + 1, m 6 + 1) 
+(777.5 — fii) G() — 771 x G(m\ + 1, m.3 — 1) + m i G{mi + 1, ms — 1) 
+7776 G(me + \)u\ + 7775 G(777 5 + 1) 77^ + mi G(tt7i + 1) u\ — uiq G(itiq + 1) u\ 

— 777 5 G(777 5 + 1) u\ + 777! (5(777! + 1) 77g + 777 6 G(t77 6 + 1) 777^ — 777 6 G(777 6 + 1) w\ 
—7775 ^(7775 + 1) w\ — 777i G(t?7i + 1) w\ . (12) 



The function G is denned as in Eq. (IA6|) . but for ease of writing only the arguments 
which contain a ±1 are shown explicitly [so that, for instance, G{m^ + 1) stands for 

(5(7771, m 2, m 3'i ^4, m 5 + 1, m &) e ^ c -} 

The general solution of these recursions is on itself of great interest, but we consider here 
only the case u\ = U2 = u% = which corresponds to Hylleraas basis set. For this we put 
mi = vtl2 = 7773 = 1 and differentiate these identities over 771,772,773 at u\ = U2 = 773 = 0. 
It leads to recursions for the h function of Eq. (jA7|) . The first subset of identities forms 
3 linear equations for h{n\ + 2),h(rri5 + l),h(rriQ + 1) (same convention for the explicitly 
written indices as for the function G) which can easily be solved. Since there are 9 equations 
for 6 unknowns, two different solutions for ^.(7774 + 1), h(m 5 + 1), h(mQ + 1) can be used to 
simplify recursion formulas. We achieve this by solving three equations 

hi (7774 + 1) - h 2 (7774 + 1) = 0, 
hi(m 5 + 1) — h,2(m 5 + 1) = 0, 

/7i(777 6 + l)-/7 2 (777 6 +l) = 0, (13) 

against the following linear combinations 

h(m 5 - 1, 777 6 + 1) - /7(777 4 - 1, 7?7 6 + 1) , 

h(m 5 + 1, m 6 — 1) — h(m A - 1, m 5 + 1) , 

h(rri4 + 1, 7775 — 1) — hijrii + 1, t?7 6 — 1) . (14) 

The obtained solutions are inserted back in recursion formulas for h(rii + 2) and himi + 1), 
and they take the form 



2w\ w\ 



rii + 



j- h{ni + 2) = [w£ (2 + 77i + n 2 + 2 773) + w\ (2 + m + 2 ti 2 + n 3 ) 



S 



w\ (2 + 2 m + n 2 + n 3 )] + w 3 Q 3 + w| Qi - w\ Q x + . . . 

2 it?^ vo^ 

l — f- h(n 2 + 2) = [vol (2 + ri! + n 2 + 2 n 3 ) - w\ (2 + m + 2 n 2 + n 3 ) + 

77-2 + 1 

(2 + 2 m + tt- 2 + n 3 )] h{) + w 2 3 Q 3 - w\ Q 2 + w\ Qi + . . . 

2 up 1 ixp 1 

3—^ /i(n 3 + 2) = [-wf (2 + n a + n 2 + 2 n 3 ) + w\ (2 + n x + 2 n 2 + n 3 ) + 

77-3 + 1 

iw? (2 + 2 77i + n 2 + 773)] h() - wlQ-i + w 2 2 Q 2 + wlQx + . . . (15) 

to? m 4 /i(m 4 + 1) = (-77x + 77 2 + tt- 3 + 2 m 4 - 1) /i() + 2 (X 3 - X 2 ) + . . . 

w 2 m 5 /i(m 5 + 1) = (m - n 2 + n 3 + 2 m 5 - 1) hQ + 2 (X l - X 3 ) + . . . 

w 3 2 m 6 /i(m 6 + 1) = (m + 77-2 - n 3 + 2 m 6 - 1) + 2(X 2 - + . . . (16) 

where by ... we denote the omitted boundary terms, which are proportional to 5 ni , and 

r7 3 + 1 77 2 + 1 

ft,(n 3 + 2, 7775 — 1) H 7 ^(w 2 + 2, m 6 — 1) , 



77 3 + 1 ' 77 2 + 1 

Q 2 = ^^h(n 2 -2,n 3 + 2) + !^^) %1 + 2,n 2 -2) + 

r7 3 + 1 77i + l 

/i(77-3 + 2, 777-4 — 1) H r fr( n l + 2, 777 6 - 1) , 



r7 3 + 1 772 + 1 

773(773-1) 723(713-1) w LO x . 

Qi = — — h(n 2 + 2, 773 - 2) H — — h{ni + 2, 773 - 2) + 

77 2 + l Til + 1 

' /i(t7 2 + 2, m 4 - 1) + — — /i(n x + 2, m 5 - 1) , (17) 



77 2 + l 72i + l 

-Xi = (n 2 - 1) 77 2 m 4 h(n 2 — 2, m 4 + 1) — (773 - 1) 77 3 m 4 ft(n 3 — 2, m 4 + 1) , 
X 2 = (n 3 - 1) 77 3 m 5 ft(n 3 — 2, m 5 + 1) — (m - 1) n\ m 5 ft,(ni — 2, m 5 + 1) , 
X3 = (771 - 1) ni m 6 h(ni — 2, m 6 + 1) — (tt- 2 - 1) tt- 2 7n 6 ft,(n 2 - 2, m 6 + 1) . (18) 

There is obviously some freedom in using Eqs. (fT5j) to simplify recursions, and our choice 
allows to separate recursions into two subshells of tt-i + n 2 + 77.3 and m 4 + m 5 + m 6 . We 
put 7774 = m-5 = 7776 = 1 in Eq. ()15|) and immediately obtain the first subset of recursions, 
Eq. (j3J). The second subset Eq. (J3J) is obtained from Eq. (|TB|) in three steps. The first step 
is setting 777-4 — m 5 — m e — 1- The second step is multiplication of recursion formulas by 
Wi Wj to remove any Wi from the denominator. The third step is differentiation with respect 
to Wi,w 2 , and w 3 , which converts the function h() into the function /(). That completes 
the derivation of recursion relations for three-electron Hylleraas integrals. 



V. SUMMARY 

We have derived a complete set of recursion relations for three-electron Hylleraas integral 
/. They are sufficiently stable for a precise numerical calculation. Since all to; are bounded 
from below by approximately \/2E, where E is the ionization energy, the denominators 
of three-electron recursions are also bounded from below. Therefore, by using extended 
precision arithmetics (octuple precision) we can safely devote 1 digit per iteration, and 
still preserve quad precision for about 30 iterations, which is the maximum we aim to use. 
3 ossible instabilities in two-electron integrals are avoided by applying inverse recursions 
loj | . Numerical results are shown in Table These recursion relations can be used for high 

TABLE I: Values of three-electron Hylleraas integral at w\ = W2 = vus = 1, [m] denotes 10 m 



n /(n, 0,0; 0,0,0) 

2.208 310 154 388 618 874 536 424[-l] 

1 2.876 820 724 517 809 274 392 190[-1] 

2 6.071 253 765 587 525 062 881 067[-l] 

3 1.801 456 579 614 247 419 513 752[0] 

4 6.949 688 537 201 117 333 162 822[0] 

5 3.316 893 553 498 521 645 367 878[1] 

6 1.892 427 697 247 010 803 401 964[2] 

7 1.258 719 915 821 483 876 136 660[3] 

8 9.575 385 319 725 442 534 735 866[3] 

9 8.206 804 555 680 135 296 239 238[4] 



/(0,0,0;n,0,0) 

2.208 310 154 388 618 874 536 424[-l] 
2.208 310 154 388 618 874 536 424[-l] 
3.658 582 716 243 175 207 969 277[-l] 
8.803 723 087 040 150 596 505 449 [-1] 
2.849 464 173 126 685 211 199 798[0] 
1.176 795 411 671 425 935 279 582[1] 
5.962 899 567 501 152 778 486 008[1] 
3.596 955 116 745 326 378 301 909[2] 
2.522 862 411 307 814 783 043 058[3] 
2.019 476 554 953 447 619 512 494[4] 



precision calculation of lithium wave function with a small computational effort. We think 
that similar recursions can be derived also for beryllium, which involves four electrons, but 
the generic integral with all powers of equal to —1 has not yet been obtained. However, 
this generic integral can be represented as a multiple sum as recently indicated by Sims and 
Hagstrom [12j and Frolov [13] and thus obtained numerically. Coming back to lithium, our 
aim is the high precision calculation of higher order relativistic and QED effects like that for 
lithium hyperfme splitting. It requires, however, also the calculation of Hylleraas integrals 
with various negative powers of and r^. 1/r" with n > 1 can be expressed in terms of 
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harmonic polylogarithms, which were introduced recently in [14], and 1/r™- can be obtained, 
we think, by generalization of the recursion relations obtained in this work. 

In summary, we have presented a simple method to derive recursion relations and solved 
them for three-electron Hylleraas integral. 
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APPENDIX A: DEFINITIONS OF FUNCTIONS AND RELATIONS 



We use the symbol 5 n which denotes the Kronecker 5 n ^. Definitions of all functions used 

— * — * 

in this work are presented below. r\j denotes |fj — fj\ and = \ki — kj\. All the rij and rrii 
are assumed to be nonnegative integer and w^Ui to be real and positive. 

d 3 r\ f d 3 r 2 f d 3 r 3 



F(n l ,n 2 ,n 3 ;n 4 ,n 5 ,n 6 ) 



An J An J An 



-wi r\—W2 T2—W3 r^—u\ r23— «2 rsi—u^ r\2 



r ni-l „n 2 -l „n 3 -l n 4 -l „n 5 -l n 6 -l 
'23 '31 '12 '1 '2 '3 



(Al) 



o? 3 ri f d 3 r 2 f d 3 r 3 



f <F r x f cT r 2 /" 
f{ni,n 2 ,n 3 ;n 4 ,n 5 ,n 6 ) = J -j— J -j— y 



«>i ri — u>2 »"2 — u>3 r3 



An 



r ni-l „n 2 -l „n 3 -l „n 4 -l „n 5 -l „n 6 -l 
'23 '31 '12 ' 1 ' ■> ' " 



(A2) 



fd 3 r 1 fd 3 r 2 f d 3 
f{*,n 2 ,n 3 ;n A ,n 5 ,n 6 ) = J -j— J -j— J 



r 3 



-wi r i — u>2 r 2 — W3 r% 



47T 7 47T 

' 1 o ' 



4 7T 5 3 (r 23 ) r™ 2 - 1 ^r 1 r " 4_1 < 5_1 r ™ 6 ™ 1 



(A3) 



/(ni,n 2 ,n 3 ;^,n 5 ,n 6 ) = 



d 3 ri f d 3 r 2 f d 3 r 3 



-Wl r\ — W2 T2 — W3 T3 



An J An J An 

„m— l ^2-1 ^,"3-1 

23 



r^rg-UTr^rOrJ-- 1 ^- 1 



(A4) 



(i 3 ri /■ d 3 r 2 



/or 7"i f 



An 



-ai n-a2 T2-Q3 ri2 r ni—l „7i2-l r n3—l 
1 1 '2 ' 12 



(A5) 



G(mi, m 2 , m 3 ; m 4 , m 5 , m 6 ) = ^ | d 3 ^ | d 3 A; 2 | rf 3 A; 3 + u \)- mi [k\ + ^)~ m2 (A6) 



{kl + ui)- m * {k 2 32 + wl)— (kf 3 + wl)— (k 2 21 + wl)'™* 



h(n 1 ,n 2 ,n 3 ;m 4 ,m 5 ,m 6 ) 



dul 



ui=0 



du7 



f-l 



\»3 



«9™ 3 



M2=0 



"3=0 



(7(1, 1, l,m 4 ,m 5 ,m 6 ) 



(A7) 



Qn 4 

f(ni,n 2 ,n 3 ;n 4 ,n 5 ,n 6 ) = (-1)™ 4 ^-^ (-l) n 



Qn 5 



l) m -g-^h(m,n 2 ,n 3 , 1,1,1) 



(A8) 
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